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Abstract 

In this brief report we compare the predictions of a recent next-to-next-to-leading order hard- 



Oh 

thermal-loop perturbation theory (HTLpt) calculation of the QCD trace anomaly to available 
lattice data. We focus on the trace anomaly scaled by T 2 in two cases: Nf = and Nf = 3. 

m 

When using the canonical value of \x = 2irT for the renormalization scale, we find that for Yang- 
t-h Mills theory (Nf = 0) agreement between HTLpt and lattice data for the T -scaled trace anomaly 

m 

begins at temperatures on the order of 8T C while when including quarks (Nf = 3) agreement 
O begins already at temperatures above 2T C . In both cases we find that at very high temperatures 



the T 2 -scaled trace anomaly increases with temperature in accordance with the predictions of 
HTLpt. 

PACS numbers: 11.15.Bt, 04.25.Nx, ll.10.Wx, 12.38.Mh 



1 



Introduction : The QCD equation of state (EoS) is of crucial importance for a better 
understanding of the matter created in the ultrarelativistic heavy- ion collisions [1J, as well 
as the candidates for dark matter in cosmology [2]. Due to the nonperturbative nature 
of QCD at low temperatures and small chemical potentials where the system is strongly 
coupled, lattice gauge theory is by far the most successful method for determining the EoS 
in this part of the phase diagram. However, due to the asymptotic freedom of QCD, lattice 
calculations should be compatible with perturbation theory at very high temperatures where 
the coupling constant becomes asymptotically weak. 

The calculation of QCD thermodynamics utilizing weakly-coupled quantum field theory 
has a long history [3H5] . The free energy is now known up to order g® log g s j5] . Unfortu- 
nately, a straightforward application of perturbation theory is of no quantitative use at the 
phenomenology-relevant intermediate coupling regime which is probed in the ongoing ultra- 
relativistic heavy-ion collision experiments at Brookhaven National Laboratory's Relativistic 
Heavy Ion Collider and the European Organization for Nuclear Research's Large Hadron 
Collider (LHC). The issue is that the weak-coupling expansion series oscillates wildly and 
shows no sign of convergence unless the temperature is astronomically high. For example, 
simply comparing the magnitude of low-order contributions to the QCD free energy with 
three quark flavors (Nf = 3), one finds that the g^ contribution is smaller than the g 2 s con- 
tribution only for g s < 0.9 or a s < 0.07 which corresponds to a temperature of T ~ 10 5 GeV 
~ 5 x 10 5 T c , with T c ~ 170 MeV being the QCD critical temperature. 

The poor convergence of the weak-coupling expansion of thermodynamic functions stems 
from the fact that at high temperature the classical solution is not described by massless 
particles, but instead massive quasiparticles due to plasma effects such as the screening 
of (chromo)electric fields and Landau damping. One way to deal with the problem is to 
use an effective field theory framework in which one treats hard modes using standard 
four-dimensional QCD and soft modes using a dimensionally reduced three-dimensional 
SU(3) plus adjoint Higgs model [SHZj but treating the soft sector non-perturbatively by 
not expanding the soft contribution in powers of the coupling constant [8]. This picture 
of treating the soft sector non-perturbatively is ubiquitous and there exist several ways of 
systematically reorganizing the perturbative series at finite temperature [9] . Such treatments 
are based on a quasiparticle picture in which one performs a loop expansion around an 
ideal gas of massive quasiparticles, rather than an ideal gas of massless particles. In this 
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brief report, we will present results for the trace anomaly in pure Yang-Mills (Nf = 0) 
and QCD with Nf = 3 obtained using the hard-thermal-loop perturbation theory (HTLpt) 
reorganization [TUT - TPS"] . We will focus on the quantity (£ — 3V)/T 2 which, in pure Yang-Mills, 
has been shown to be approximately constant for temperatures which are close to the phase 
transition temperature. We will show that, when using the trace anomaly scaled by T 2 , 
(£ — 3V)/T 2 , as a criterion, pure Yang-Mills exhibits a transition to perturbative behavior 
at temperatures on the order of 8 T c while for Nf = 3 one sees reasonable agreement with the 
next-to-next-to-leading order (NNLO) HTLpt result already at temperatures on the order 
of 2-3T c . 

Background : HTLpt is a reorganization of the perturbative series for thermal gauge 
theory that is a gauge-invariant generalization of screened perturbation theory (SPT) applied 
to scalar 4 theory [13] . In scalar 4 theory the basic idea of SPT is to add and subtract a 
thermal mass term from the bare Lagrangian: including the added piece in the free part of the 
Lagrangian; while treating the subtracted piece as an interaction term on the same footing 
as the quartic term. In gauge theories, however, simply adding and subtracting a local mass 
term violates gauge invariance [2]. Instead, one adds and subtracts the hard-thermal-loop 
(HTL) effective action, which dresses the propagators and vertices self-consistently so that 
the reorganization is manifestly gauge invariant [15J. HTLpt thermodynamic calculations 
for non-Abelian theories have been recently carried out to three loops, aka NNLO [Hi [12] 
and improved convergence comparing to the weak-coupling expansion is found. 

At temperatures just above the phase transition, 1.1 T c < T < 4T C , in pure Yang-Mills, 
lattice results find that the trace anomaly scaled by T 2 is approximately constant. It has been 
suggested that this behavior is due to the influence of power corrections of order T 2 [16] which 
are nonperturbative in nature and might be related to confinement. Phenomenological fits of 
lattice data which include such power corrections show that the agreement with lattice data 
is improved [17] . Alternatively, the power corrections are obtained from AdS/QCD models 
which break conformal invariance [T8TI2T] . The potential impact of the power corrections on 
relativistic hydrodynamics relevant to heavy ion experiments at LHC is explored in Ref . [22J . 
However, it is important to mention that current AdS/QCD models have no firm theoretical 
reason for including modifications which induce power law corrections besides their seeming 
phenomenological importance. In addition, the models are largely unconstrained in the 
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NNLO HTLpt - 1 loop running, /i = nT 

NNLO HTLpt - 1 loop running, yu = 2nT 

NNLO HTLpt - 1 loop running, yu = 4nT 
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FIG. 1. NNLO HTLpt result for the pure Yang-Mills (Nf = 0) trace anomaly over temperature 
squared compared with lattice data from the Wuppertal-Budapest group [21] at high temperatures. 
For all curves we have used the one- loop running of the strong coupling as specified in the text. 
Different HTLpt curves correspond to different choices of the renormalization scale /i. 

magnitude of these corrections and the parameters in the underlying models must be fit to 
lattice data. 

In the last decade there have been intensive studies of the EoS on the lattice for both 
pure Yang-Mills and full QCD with dynamical quarks. These studies generally measure the 
trace anomaly scaled by T 4 , {£ — 3V)/T 4 , and use this to determine the energy density, 
pressure, and entropy density of the system. In the remainder of this brief report, we 
compare the NNLO HTLpt trace anomaly with lattice data from the Wuppertal-Budapest 
[21125], HotQCD [261 [27], RBC-Bielefeld [281 [29], and WHOT-QCD [301 HQ collaborations. 

Results : In Fig. [I] we plot the trace anomaly scaled by T 2 T C 2 for pure Yang- Mills theory 
and compare this with recent lattice results obtained by the Wuppertal-Budapest collabo- 
ration [21]. This figure shows the behavior of A = (S — 37 , )/(T 2 T c 2 ) at high temperatures. 
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FIG. 2. NNLO HTLpt result for the pure Yang-Mills (Nf = 0) trace anomaly over tempera- 
ture squared compared with lattice data from the Wuppertal-Budapest [21], Bielefeld [23], and 
WHOT-QCD [301 ED] collaborations. For the NNLO HTLpt result with three-loop running, the 
renormalization scale was taken to be \i = 2nT. For the NNLO HLTpt result with one-loop 
running, we show a line corresponding to fx = 2nT and a band showing the range ttT < fj, < AttT. 

For the NNLO HTLpt result we use the perturbative one-loop running of the strong cou- 
pling which is consistent with the counterterms that one obtains at NNLO. We show three 
different curves corresponding to the three different values for the renormalization scale, 
/i = {itT, 2ttT, AttT}. The central value of = 2ttT is the scale one expects to emerge from 
an equilibrium calculation since the lowest possible non-vanishing bosonic Matsubara mode 
sets the scale; however, for completeness sake we show the variation of /i by a factor of two 
in either direction, which should be understood as a dramatic overestimate of the theoretical 
uncertainty. As can be seen from Fig. [TJ at low temperatures NNLO HTLpt underpredicts 
the trace anomaly for pure Yang-Mills and only starts to agree at temperatures on the order 
of 8 T c for fi = 2nT. The agreement with lattice results begins at much higher temperatures 

5 



for the other extreme values of /i. At asymptotically high temperatures we see excellent 
agreement between the NNLO HTLpt results and the lattice results. We note in this con- 
text that if one allows for a fit to the unknown perturbative g 6 contribution to the pressure, 
then the EQCD approach can equally-well describe the pure Yang-Mills lattice data down 
to temperatures on the order of 8T C [6j [7]. In the EQCD case, once again, the key is to not 
expand perturbatively in the Debye mass, but to instead treat contributions from this scale 
non-perturbatively. 

The behavior of the trace anomaly at lower temperatures is presented in Fig. [2] along with 
lattice data from the Bielefeld [23J and WHOT-QCD [301131] collaborations. The solid black 
line is the NNLO HTLpt result obtained using a one-loop running of the QCD coupling 
and the dashed black line is the HTLpt result obtained using a three-loop running of the 
QCD coupling. In the case of the three-loop running we use the lattice determination of 
T c /A^g = 1.26 to fix the scale [24]. For comparison between the one- and three-loop results 
we require that both give the same value for the strong coupling when the renormalization 
scale /i = 5 GeV. Numerically, one finds a s (5 GeV) = 0.140553. We compare one- and 
three-loop running because, formally speaking, one only obtains counterterms consistent 
with one-loop running from the NNLO HTLpt calculation; however, at temperatures close 
to the critical temperature the three-loop running gives important corrections to the scale 
dependence of the running coupling [32] • We will use the difference between one- and three- 
loop running to gauge the theoretical uncertainty of the NNLO HTLpt results. For both 
the one- and three-loop running the lines assume \i = 1nT to fix the renormalization scale 
as a function of the temperature. For the one-loop running case we also show a band 
corresponding to varying the renormalization scale in the range 7rT < \i < AttT. We note 
that this value is not fit in any way to the lattice data, but instead is the predicted value of 
the appropriate renormalization scale to use. Better agreement at high temperatures can be 
obtained by fitting \i to lattice data, in which case one finds that a value of fi/ (2nT) = 1.75 
is preferred in the case of pure Yang-Mills [24] . Here we present the results at face value 
without fitting. 

The most remarkable feature of the lattice data for A is that it is essentially constant in 
the temperature range T ~ 1.1 — 4T C . At temperatures above approximately 4T C the latest 
Wuppertal-Budapest [2H [25] results show an upward trend in accordance with perturbative 
predictions. The WHOT-QCD results also exhibit an upward trend, however, it starts at 
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FIG. 3. NNLO HTLpt result for the Nj = 3 QCD trace anomaly over temperature squared 
compared with lattice data from the Wuppertal-Budapest |25j, hotQCD [26j[27], and RBC-Bielefeld 
|28|, 129] collaborations. For both HTLpt curves the renormalization scale was taken to be fj, = 2nT. 



much lower temperatures. This discrepancy could be due to their fixed scale approach not 
having sufficiently large N T at high temperatures as noted in their paper [30] . 

In Fig. H we present the N f = 3 QCD NNLO HTLpt prediction for (£ - 3V) /T 2 and 
compare to lattice results available from the Wuppertal-Budapest [25] and hotQCD col- 
laborations [26l E7] . The lattice data from the Wuppertal-Budapest collaboration uses the 
stout action and have been continuum estimated by averaging the trace anomaly measured 
using their two smallest lattice spacings corresponding to N T = 8 and N T = 10. The lattice 
data from the hotQCD collaboration are their iV T = 8 results using the asqtad, p4, and 
HISQ actions which have not been continuum extrapolated [26, 27J. The lattice data from 
the RBC-Bielefeld collaboration is iV r = 6 and have also not been continuum extrapolated 

[2HIE9]. 

As before, we present HTLpt results using both one- and three-loop running of the 
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FIG. 4. NNLO HTLpt result for the Nj = 3 QCD trace anomaly over temperature squared 
compared with lattice data from the Wuppertal-Budapest |25j, hotQCD [261 [27], an( l RBC-Bielefeld 
|28|, [29] collaborations. For the HTLpt curve we use the one-loop running of the strong coupling 
constant. The red band is the variation of the HTLpt result with the assumed renormalization 
group scale (ttT < fi < AttT) and the line corresponds to the value /U = 2irT. 



strong coupling. In this case we have required a s (5 GeV) = 0.2034 in accordance with 
recent high precision lattice measurements of the running coupling constant [33]. As can 
be seen from Fig. [3] for T > 400 MeV we find reasonable agreement between the NNLO 
HTLpt predictions and available lattice results. Below this temperature the Wuppertal- 
Budapest and hotQCD results do not seem to agree. Therefore, it is difficult to draw 
conclusions about the efficacy of the HTLpt approach; however, naively one expects that 
for temperatures less than twice the critical temperature that non-perturbative corrections 
should become increasingly important. 

In Fig.Qwe again present the N f = 3 QCD NNLO HTLpt prediction for (£-3V)/T 2 and 
compare to lattice results. However, in this case we show only the result of assuming one-loop 
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running for the strong coupling constant which is consistent with the counterterms necessary 
to renormalize HTLpt at NNLO. The red band results from varying the renormalization 
scale /i by a factor of two around /i = 2-kT which is the theoretically favored value. As 
can be seen from this figure, it would seem at first glance that there is still a significant 
theoretical uncertainty related to the choice of the renormalization scale /i; however, fi is 
not a completely free parameter since its value should be related to the lowest possible 
momentum exchange, which for a plasma in equilibrium is set by the lowest non- vanishing 
bosonic Matsubara mode. In this sense, the band shown dramatically overestimates the 
theoretical error of the final result; however, we include it here for completeness. 

Discussion and Conclusion : The purpose of this brief report is to present an easily 
accessible comparison of predictions of the T 2 -scaled trace anomaly with recent NNLO cal- 
culations using HTLpt. The details of the HTLpt framework and the full NNLO calculation 
itself can be found in Refs. [TTJ [TJ]. We compared the N c = 3 NNLO HTLpt predictions 
with the QCD T 2 -scaled trace anomaly in two specific cases, namely Nf = and Nf = 3, 
however, the result for general N c and Nf can be found in Refs. (TTJ [12]. 

We find that by including quarks in the calculation, agreement with lattice data is greatly 
improved as compared to the NNLO results of pure-glue QCD. Fermions are perturbative in 
the sense that they decouple in the dimensional-reduction step of effective field theory, so we 
expect that including contributions from quarks gives at least as good agreement with the 
lattice calculations as compared to the pure-glue case. However, the reasons for the large 
improvement between the HTLpt predictions and the lattice calculations is not clear to us. 

The T 2 -scaling of the trace anomaly has received much attention recently and comparisons 
with lattice data provides an extremely sensitive test of the agreement between theory and 
numerical results. Using \x = 2ttT we find that for Yang-Mills theory (Nf = 0) agreement 
between HTLpt and lattice data for the trace anomaly begins at temperatures on the order 
of 8 T c while when including quarks (Nf = 3) the agreement begins already at temperatures 
above 2 T c . In both cases we find that at very high temperatures the T 2 -scaled trace anomaly 
increases with temperature in accordance with the predictions of HTLpt. 

In closing we emphasize again that the T 2 -scaled trace anomaly is a very sensitive ob- 
servable. If one were to compare instead the pressure, energy density, and entropy density 
one finds that, in both cases presented here, HTLpt is consistent with lattice data down 
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to temperatures on the order of 2 — 3 T c (TTJ [12] . Looking forward, we point out that the 
HTLpt framework is formulated in Minkowski space and is therefore also applicable to the 
calculation of real-time properties of the quark gluon plasma at LHC temperatures. 
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